Elasticity of arrested short-ranged attractive colloids: 
homogeneous and heterogeneous glasses 
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We evaluate the elasticity of arrested short-ranged attractive colloids by combining an analytically 
solvable elastic model with a hierarchical arrest scheme into a new approach, which allows to dis- 
criminate the microscopic (primary particle-level) from the mesoscopic (cluster-level) contribution 
to the macroscopic shear modulus. The results quantitatively predict experimental data in a wide 
range of volume fractions and indicate in which cases the relevant contribution is due to mesoscopic 
structures. On this basis we propose that different arrested states of short-ranged attractive colloids 
can be meaningfully distinguished as homogeneous or heterogeneous colloidal glasses in terms of the 
length-scale which controls their elastic behavior. 
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Solutions of short-ranged attractive colloidal particles 
are the object of intense study due to their technologi- 
cal applications (proteins, paints etc.), as the constituent 
blocks of nanomatcrials, as well as model systems to in- 
vestigate phase behavior and dynamical arrest of con- 
densed matter A landscape of phases has been ob- 
served upon varying the volume fraction <f> or the inter- 
action parameters [2j, and, as a matter of fact, extended 
regions of the phase diagram are still poorly understood. 
In very dense suspensions (</> > 0.5) the arrested states 
are spatially homogeneous (i.e. the typical linear size of 
structural heterogeneity is smaller than the particle di- 
ameter Rq). Particles are immobilized within the range of 
attraction, giving rise to bonds that are persistent under 
strain in the linear regime, and the high density leads 
to (attractive) glassy states @. For 0.2 < < 0.5, 
the situation is complicated: arrested states can only 
occur thanks to pronounced structural heterogeneities, 
typically on length scales larger than i?o- Therefore, 
they are more related to gelation rather than to the 
caging typical of crowded random media [f|. However, 
such arrested states are also hardly classifiable as clas- 
sic network gels in view of the different morphology and 
stress-bearing mechanisms. For them, different micro- 
scopic phenomena should be considered and a new theo- 
retical framework, able to account for the strong spatial 
heterogeneities and currently still missing, would be de- 
sirable. A crucial point, mostly neglected in recent stud- 
ies, is that arrested states occurring at different volume 
fractions and attraction strengths do display dramati- 
cally diverse mechanical and rheological properties Q. 
That is why, their characterization is of true interest in 
technological applications and material design. 

In this Letter we propose a new, more down to earth, 
approach to characterize arrested short-ranged attractive 
colloids, based on their mechanical response. Being the 
problem extremely hard to tackle, we rely on a simpli- 
fied picture: We combine an analytically solvable elastic 



model with a hierarchical arrest scheme and provide a 
first attempt to discriminate the microscopic (primary 
particle-level) from the mesoscopic (cluster-level) contri- 
bution to the macroscopic elasticity. The predictions are 
given in terms of the shear modulus, a quantity that can 
be measured in a rheometer, and are directly compared 
with experimental results. In spite of the crude, but phys- 
ically grounded assumptions, we are able to show, for 
the first time, that where structural heterogeneities as 
clusters are present, the measured macroscopic elasticity 
differs by orders of magnitude from the one of the ho- 
mogeneous glass because it is strongly dominated by the 
inter-cluster contribution. Our approach offers therefore 
a new insight into the Physics of arrested attractive col- 
loidal suspensions, especially in the range 0.2 < </> < 0.5, 
hitherto poorly understood. We propose that homoge- 
neous and heterogeneous arrested states can be distin- 
guished in terms of the length scale i?o, which controls 
the macroscopic elastic modulus. 

Model and predictions. We study a suspension of col- 
loidal particles interacting via a short-range attraction, 
as the one typically induced by depletion using non ad- 
sorbing polymer well above percolation, so that a 
finite shear modulus is always detectable. As a first 
step towards rationalizing the mechanical response of ar- 
rested, dense attractive colloids, and in the fundamen- 
tal lack, at present, of alternative approaches, we use a 
simplified, although physically grounded, scenario. We 
consider that, by reducing the volume fraction from the 
very dense homogeneous glass case, the attractive inter- 
actions will start to produce small aggregates, with the 
aggregation being soon arrested by the glassy environ- 
ment. This mechanism leads to compact aggregates of 
relatively small linear size 0, I n order to distin- 
guish the single-particle contribution to elasticity from 
the mesoscopic one, we therefore consider coarse-grained 
entities, i.e. renormalized particles. Their effective inter- 
actions are directly responsible for the macroscopic prop- 
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erties of the system. This scheme can be associated to a 
double crgodicity-breaking scenario where the arrest oc- 
curs in form of a cluster-glass transition sine 
qua non is that further coalescence of the clusters, lead- 
ing to phase separation, is prevented on the observation 
time scale ■ We distinguish between the macroscopic 
elasticity of the system G, the intra-cluster elasticity G g 
resulting from mutual interactions between primary par- 
ticles, and the inter-cluster elasticity G c resulting from 
mutual interactions between clusters. 

Let us first consider the case of a homogeneous at- 
tractive glass in the low-temperature and high-density 
((f) > 0.5) region of the phase diagram. Here 
we focus on the first linear regime reported in rhe- 
ological measurements |l2|, which is due to bond- 



breaking and, for strong attractions, extends nearly up 
to strains of a few percent. Hence, we follow the 
Cauchy-Born approach developed in [13j for amorphous 
solids and obtain the elasticity from a free energy ex- 
pansion around a stressed (quenched) reference state 
{R} where all particles are labeled. The expansion 
reads SF = F({r}) - F({R}) ~ [dF/dr l3 ]Sr l3 + 
\ Ylij,ki [d 2 F/dr i:j dr k i] SrijSr k i, where the sums run over 
all interacting pairs of particles (bonds) and the deriva- 
tives are evaluated at relative distances {Ry} in the 
quenched reference state. Due to the purely internal 
nature of the stresses, the first term in the r.h.s. does 
not contribute to the macroscopic elasticity 13|, [l4[ . In 
the case of central forces, considering a pair-interaction 
with a deep minimum e 3> fcsT, the bond stiffness is de- 
fined by k — d 2 F/drfj, and the Born-Huang (BH) terms 

become 5F = 5^5^y(^ r »j) 2 - ^ we express such devi- 
ations in terms of a smooth displacement field Uy [HI, 

with central forces only the longitudinal component ujj is 
non-zero so that, upon neglecting higher order terms, the 
deformation free energy reads SF ~ |*sEy(«{,-) 2 - Un- 
der the assumption that all pairs of particles are much 
localized near the minimum of the potential well and 
that displacements are affine, we obtain the continuum 
limit of the microscopic u\- in terms of the macroscopic 
strain tensor e, (ufj) = Tr[(Ry <& Ry) ■ e/Rij], with 

(4-<4» a < <4> 2 ' toobtain(5F - ^Ei^KRy® 
Rij) • e/Rij]} 2 . For the quenched reference configuration 
{R}, in the case of uncorrelated disorder, the summation 
over bonds can be replaced by the total number of con- 
tacts, which introduces the average coordination number 
z((j)) in the final expression of the deformation free en- 
ergy. We then derive the off-diagonal components of the 
stress tensor 1151 
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where (•) f2 = (Air) 1 / •dCl. Hence the shear modulus 
of an arrested phase of volume fraction <fi = P^Ro/ft, 



with mean particle diameter ~ Rq, is given by G g = 
(2/5)ir- 1 (f>z((t>)KR 2 l - d . The evaluation of the coordina- 
tion number z(<p) in dense glassy systems is still contro- 
versial. Since at high densities also attractive systems 
are likely to be dominated by excluded-volume repul- 
sion [l7[, we estimate the number of short-range con- 
tacts involved in the mechanical bonds from z(<f>) of a 
hard sphere glass. That is, we integrate the radial dis- 
tribution function (rdf) of the dense (hard-sphere) liquid 
precursor, g(r), within a shell of width ft . A reasonable 
choice for repulsive systems is ft ~ 0.03 which allows to 
recover z ~ 6 at random close packing [l6[. We thus 
obtain the shear modulus as a function of the rdf: 



G g ~ (48/5) 7 r- 1 Ki?2- d / (l + lf g (l- <f>)dl (2) 
Jo 

where I = (r - R )/R and ft ~ 1/30. For g(r) near con- 
tact (I < 0.1) we use liquid theory valid in the dense hard- 
sphere fluid [l9| and we calculate k using the Asakura- 
Oosawa (AO) potential [2(j. Our predictions have been 
compared to the experimental data of Ref. [l2| for an at- 
tractive colloidal glass (4> ~ 0.6). The shear modulus was 
measured as a function of the attraction strength, i.e. the 
reduced polymer concentration c p /c*. For strong enough 
attraction (and fully elastic response), the affine approx- 
imation is more realistic and our prediction gives an ac- 
curate estimate of the measured shear modulus [HI . We 
have checked that the agreement does not significantly 
change upon vary ft . 

For a more stringent test of the model where the bond- 
stiffness and the local structure of the glass are well- 
defined, we performed Molecular Dynamics simulations 
of deeply quenched glasses obtained from different mod- 
els for supercooled liquids. We subjected the glass to 
shear strain and extracted stress-strain curves [18|. The 
agreement obtained in the linear regime is remarkable, in- 
dicating that the model is quantitatively predictive and 
that the nonaffine rearrangements are actually small in 
the strong-attraction or deep-quench limit in agreement 
with our assumptions [l8| . 

For the more dilute regime, we now consider local ag- 
gregation of the colloidal particles to form beads changed 
by small variation of ft. (clusters) which, in turn, 
arrest due to cither caging or residual attraction, as 
in the double-ergodicity breaking scenario mentioned 
above 

El EH- The clusters are viewed as compact 
(spherical or quasi-spherical) renormalized particles of di- 
ameter Rq, whose effective volume fraction may be iden- 
tified with the one determined by the spheres enclosing 
them (i.e. significantly larger than <f). If the cluster lin- 
ear size is larger than the particle diameter by a factor 
say less than 10, each contact between clusters is likely 
to reduce to a single colloid-colloid bond. Upon neglect- 
ing: 1) the breakup probability within the cluster, and 2) 
the effect of long-range repulsion, the effective interaction 
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FIG. 1: (color online) Shear modulus as a function of <f> c for 
the experimental data from [?J and for the model prediction 
from Eqs.Q and @- RqIRo — 5.5 and polymer-to-colloid 
size ratio £ = 0.078 from [7J (upper x-axis: <j> for c p /c* — 0.3). 
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FIG. 2: (color online) The shear modulus as a function of 
Cp/Cp from Ref. 0] (symbols), with polymer-to-colloid size ra- 
tio £ = 0.08 and Rq/Rq from experiments. Lines: model 
prediction from Eq.(j2j (dashed) and ((3} (full). 



between clusters obviously reduces to the bare colloid- 
colloid interaction, in agreement with [icj |. Further, the 
mean coordination will change to z(<f> c ) (where now <p c is 
the cluster volume fraction), but its form can be still de- 
termined as in Eq.© if <j) c is in the dense glassy regime 
dominated by mutual impenetrability. Hence, for the 
modulus of the cluster-glass we can write: 

G c ~ (2/5)n- 1 ^ c z(^ c )kR Q 2 ' d (3) 

with k — k under the assumption of small clusters. Eq. 
([3]) gives the elastic modulus of the material, where the 
macroscopic elasticity is dominated by the mesoscopic 
level. We test our scheme using the extensive experi- 
mental data of Ref. Q for a system of colloidal silica 
particles with polystyrene as depletant in organic solvent 
(decalin) , where the ratio of the polymer gyration radius 
to R a is the same as in the experiments of [l2j ■ For two 
values of c p /c*, they vary <f> in the range 0.2 — 0.4. The 
system is a dense gel of compact clusters (</> ~ 0.5 inside 
the cluster), whose diameter Rq/Rq ~ 5 is determined 
from the Debyc-Bueche plot. For a fixed c p /c*, they find 
Ro/Ro and the volume fraction inside the clusters not 
to significantly change with </>, i.e. <j) c increasing linearly 
with <f>. This indicates that in this case our assumptions 
on the structure are reasonable. The assumption that the 
spheres enclosing the clusters are densely packed yields 
4> c ~ 0.64 at the largest value of </>, allowing for a physi- 
cally reasonable estimate of <f> c in the range of <fi consid- 
ered. Using the experimental measurement of Rq/Rq we 
can therefore calculate G c from Eq. (J3]) . In FigQ]we plot 
the values of the shear modulus measured at 1Hz in Q 
for two different attraction strengths, together with our 
predictions of G g from Eq.® and of G c from Eq.©, 



as function of <fr and 4> c . The excellent agreement between 
our estimate and experiments in FigfT] clearly indicates 
that, differently from the case previously considered, now 
is G c that dominates the macroscopic elasticity of the sys- 
tem. We have checked that the uncertainty on <fi c does 
not significantly affect this result. For lower attractions, 
the system will gradually cross over towards the repulsive 
case. 

We now consider the results recently reported in 
for a different arrested AO system in the heterogeneous 
glass regime, studied at a fixed <f> = 0.4 and upon varying 
Cp/c* ~ 2. In the experiments, they used light-scattering 
and microscopy to estimate the length-scale of structural 
heterogeneities up to c p /c* ~ 2. We have used this ex- 
perimental input to calculate the prediction of our model 
for this system. Again for densely packed clusters (i.e. 
with volume fraction in the range 0.58-0.64) we estimate 
4> c ~ 0.6 (and check that variation around this value 
does not significantly change our outcomes). In FigfJ] 
we compare G c to the experimental measurements: G c 
(solid line) shows a remarkable agreement with G exp and 
turns out again to dominate the elastic response (see in 
comparison G g , dashed line). Also here the comparison 
with experiments has been done at low frequencies and 
for fully clastic response. 

The results just discussed lead us to sketch in FigJ3] a 
new qualitative phase diagram for the arrested states of 
short-range attractive colloidal suspensions. As a func- 
tion of <f> and IsbT/ e, where e is the depth of the attractive 
well, we locate the arrested states in the region where the 
system displays a non-zero measurable elastic modulus G 
(the continuum line divides fluid from solid states as dis- 
cussed e.g. in [23}). At high enough <f> and large attrac- 
tions, such state (homogeneous glass) is characterized by 
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FIG. 3: (color online) Arrested states characterized in terms 
of the typical length scale Ro / Ro of the mechanically relevant 
structural heterogeneities. Color gradient from R /Ro 2> 1 
(light) to Ro/ Ro — 1 (dark). Fluid-solid transition (solid line) 
and spinodal (dotted line) also indicated. 

an elastic modulus dominated by the inter-particle elas- 
ticity (Ro/Rq = 1). Upon lowering (j), aggregation starts 
to produce mesoscopic structural heterogeneities and the 
mechanical response of the system crosses over towards 
a regime dominated by the inter-cluster elasticity (het- 
erogeneous glass) (Ro/Rq 3> 1). The two regimes will 
be distinguished by a significant variation of the elastic 
modulus of the material (FigfTJ. Our diagram suggests 
for the first time a distinction between arrested states of 
attractive colloidal suspensions in terms of a well defined, 
directly measurable quantity, the elastic shear modulus. 
When the structure of the material is known in some de- 
tail, such distinction can provide information on which 
part of the structure is relevant to the elastic behavior. 
Upon further lowering the volume fraction ((/> < 0.2), we 
expect the elastic response to be dominated by a different 
length scale, associated to the weakly connected network- 
like mesoscopic structure and strongly dependent on <fi 
Due to lower connectivity, nonaffinity is likely to be 
more pronounced and lead to a sensibly lower shear mod- 
ulus (Fig[T]). At very low cj> arrested states may be due 
to effective directional interactions arising at mesoscopic 
length scales, leading to open network structures [241 ] . 
where buckling strongly affects the macroscopic elastic- 
ity. 

Conclusions. We have proposed that dramatic dif- 
ferences observed at different volume fractions in the 
mechanical response of arrested, dense attractive col- 
loids originate from the presence of structural hetero- 
geneities. We have presented a first attempt to estimate 
the shear modulus of glassy systems where structural het- 
erogeneities as clusters are present. Being the problem 
extremely hard to tackle, we have to rely on a rather 
simplified picture. In spite of this, we have been able 
to show in two different experimental systems that the 
measured macroscopic elasticity is strongly dominated by 
the inter-cluster contribution. This first novel approach 



can pave the way to a meaningful distinction between 
different glassy states, on the basis of the length-scale 
Ro dominating their elastic response. Ro varies with vol- 
ume fraction and attraction strength due to the arising of 
structural heterogeneities. However, whereas the charac- 
terization of structural heterogeneities requires detailed 
structural information and is therefore often elusive, the 
variation of Ro is unambiguously signalled by significant 
variation of the elastic modulus. Therefore we propose 
that homogeneous and heterogeneous arrested states can 
be distinguished in terms of the length scale -Ro, which 
controls the macroscopic elastic modulus. 
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